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ABSTRACT 


The  power  transmission  and  reflection  coefficients  ere  derived  for  electro¬ 
magnetic  waves  propagating  through  a  plasma  slab  having  a  parabolic  electron 
distribution.  The  analysis  considers  transverse  electromagnetic  waves  (TE 
Mode)  impinging  on  a  plasma  slab  at  arbitrary  angles  of  incidence.  The 
solutions  are  in  terms  of  complex  hypergeometric  functions  and  their  deriva¬ 
tives.  Nisnerical  results  for  the  transmission  and  reflection  coefficients 
are  plotted  as  functions  of  peak  plasma  frequency,  peak  collision  frequency, 
signal  frequency,  slab  thickness,  and  angle  of  incidence.  The  results  of 
this  study  can  be  applied  to  transmission  of  electromagnetic  energy  through 
laboratory  plasmas  that  are  bounded  by  walls.  Numerical  resuJts  are  in 
agreement  with  experimental  results  for  a  rectangular  glow  discharge  plasma. 
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SECTION  I 
INTRODUCTION 


The  propagation  characteristics  of  electromagnetic  vaves  traversing  a 
bounded  homogeneous  plasma  at  various  angles  of  incidences  are  well  known  and 
relatively  simple  to  calculate.  However,  in  most  cases  (laboratory  plasmas, 
plasma  sheaths  surrounding  reentry  vehicles,  etc.),  the  electron  density 
distribution  of  the  plasma  must  be  considered.  Many  laboratory  plasmas,  created 
by  such  devices  as  glow  discharges,  R>F  generators,  and  arc  discharges,  have 
electron  densities  that  can  be  approximated  by  a  parabolic  distribution 
(references  2,  8,  9).  Results  of  recent  experiments  using  a  plasma  hollow- 
cathode  glow  discharge  tube  indicate  that  a  parabolic  electron  density 
distribution  is  a  good  approximation  to  the  plasma  profile  (reference  4). 

In  this  present  study  an  investigation  is  made  of  the  transmission 
characteristics  of  transverse  electromagnetic  waves  propagating  through  a 
plasma  slab  having  a  parabolic  electron  distribution.  Here,  the  normal 
incidence  parabolic  case  considered  by  Bower  (reference  3)  is  extended  to 
include  arbitrary  angles  of  Incidence.  Several  errors  were  detected  in  Bower’s 
study  which  made  his  results  for  the  power  transmission  and  reflection 
coefficients  incorrect. 
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SECTION  II 
PLASMA  PROPERTIES 

In  studying  the  interaction  between  electromagnetic  waves  and  plasmas, 
only  the  effects  of  the  free  electrons  In  the  plasma  need  be  considered. 

The  mass  of  the  electrons  is  much  smaller  than  that  of  the  positive  Ions; 
therefore,  the  electrons  oscillate  about  their  equlllbrlira  position  at  a  much 
higher  frequency  than  the  ions.  When  this  frequency  Is  In  the  neighborhood 
of  the  electromagnetic  signal  frequency,  the  electrons  interact  with  the 
electromagnetic  wave  and  absorb  some  of  the  propagating  energy.  Therefore, 
the  electrical  propertl-^s  of  the  plasma,  as  "seen"  by  the  electromagnetic 
wave,  are  primarily  a  function  of  the  free  electrons  in  the  plasma. 

The  electrical  properties  of  a  plasma  slab  can  be  described  by  a  dielectric 
constant  and  an  electrical  conductivity.  These  are  functions  of  the  plasma 
frequency.  The  free  electrons  In  a  plasma  are  in  a  state  of  constant  agitation 
and  each  electron  oscillates  about  its  equilibrium  position.  The  frequency 
of  the  electron  oscillation  is  called  the  plasma  frequency  and  is  defined  as 

1  /2 

rNg(y)e^'| 

w  (y)  *  1 - I  radians /second  (1) 

P  I  "o“  J 


where 

»  electron  density  Electrons /cm 
*  permittivity  of  free  space 
m  ■  mass  of  an  electron 
e  ■  charge  of  iin  electron 

When  the  values  for  the  various  physical  constants  are  Inserted  Into  Equation 
(1),  the  plasma  frequency  Is  given  by 

u)p(y)  *  ^3.18  X  10^  radians/second  (2) 
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As  the  electrons  oscillate,  many  of  them  will  undergo  collisions  with  neutral 
particles.  The  average  number  of  collisions  per  unit  time  is  defined  as  the 
collision  frequency  The  collision  frequency  is  a  function  of  the 

electron  collision  cross  sections  of  the  various  constituents  of  the  plasma 
and  their  densities.  The  collision  frequency  is  given  by 

V  ■  I  N  Q  collisions /second  (3) 

c  ill 


where 

V  ■  average  electron  velocity 

■  density  of  each  i  constituent 

■  electron  collision  cross  section  of  the  i  constituent  (m^) 

In  this  work,  the  collision  frequency  is  assumed  to  be  constant  throughout  the 
plasma  slab.  Albini  and  John  (reference  1)  have  shown  that  the  assumption  of 
a  constant  collision  frequency  introduces  very  small  errors  when  the  electron 
density  variation  is  large  compared  to  the  variation  in  collision  frequency. 

For  a  reentry  plasma  sheath,  the  electron  density  can  vary  from  zero  at  the 
boundaries  to  10^^  electrons/cm^  and  higher  in  the  boundary  layer.  The  collision 
frequency,  however,  v^ill  usually  vary  only  by  a  factor  of  100.  In  a  hollow 
cathode  glow  discharge,  the  electron  density  can  vary  by  a  factor  of  10^^ 
across  the  plasma  with  only  a  small  variation  in  the  collision  frequency 
distribution  (reference  2). 

The  conductivity  and  permittivity  of  a  plasma  can  be  written  in  terms  of 
the  plasma  frequency,  signal  frequency,  and  collision  frequency  as,  respectively. 


a(y)  -  -J 


a)^(y)e 
p  o 

(“-jVc) 


(4) 


.•<y)  -  ^ 


(5) 


Equations  (4)  and  (5)  can  be  separated  into  real  and  imaginary  parts  to  give 


a(y) 


0  c  p  1  ^  o  _P_ 

0)2  4-  V  2  J  0)2  4  V  2 


(4a) 
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SECTION  III 

ELECTROMA^’ETIC  WAVES  IN  A  ILASMA 

Maxwell's  eoufttiofts  for  a  time-harmonic  electromagnetic  wave  of  radian 
frequency  w  in  an  isotropic  plasma  are  given  by 


V*E  - 

0 

(6) 

V«H  « 

0 

(7) 

VxE  - 

-  j 

(8) 

- 

(o  +  jwe)E 

(9) 

The  wave  equation  can  be  derived  by  taking  the  curl  of  Equation  (8),  Inserting 
the  result  into  Equation  (9),  and  using  Equation  (6).  The  wave  equation  for 
the  electric  field  in  a  plasma  is 


V^E  +  -j  E  »  0 

Equation  (10)  can  be  written  as 


we.  J 


(10) 


V^E  +  — -  E  -  0 

o  e 

0 


(11) 


where 


k  *  w/c 
o 


e:*(y) 

e 


1  -1 
•'  we 

0 


However,  the  index  of  refraction  of  a  plasma  mediim  is  equal  to  the  square  root 
of  its  dielectric  constant  (reference  7).  Therefore,  Equation  (11)  becomes 


V^S  +  k2  n2(y)  E  =  0 


(12) 
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where 


n(y) 


/ 


!£iS 


«  Index  of  refraction 


Aeame  thvtt  an  electrooagnetlc  wave  impinges  upon  an  inhoicogeneous  plasma 
slab  as  shown  in  figure  1.  Only  the  TE  mode,  for  which  the  electric  field  is 
perpendicular  to  the  plane  of  incidence,  will  be  considered.  The  angle  6^  at 
which  the  transmitted  wave  exits  the  plasma  and  the  angle  6^  of  the  reflected 
wave  will  be  the  same  as  the  angle  6^  of  the  Incident  wave.  However,  the  exit 
point  will  not  be  in  line  with  the  incident  signal  since  refraction  occurs  in 
the  plasma  (reference  7). 


Figure  1.  Geometry  of  Wave-Plasma  Interaction 
The  wave  equation  can  now  be  written  as 


32e  (y.z)  a^E  (y.z) 

— 3^2—  +  — 3ir-  +  k^'(y)E^(y.z)  -  o 


(13) 


The  Maxwell  equations  relating  the  curl  of  E  with  H,  Equation  (8) ,  are 

dE 


and 


ja)P  H  »  +  ^ 

0  z  8y 


3E 

Jcoij  H  *  -  T— ^ 
0  y  9z 


(14) 


(15) 


6 


ApyL-TR-66-70 


A  solution  to  Equation  (13)  by  the  method  of  separation  of  variables  (see 
Appendix  I)  Is 

jk  zslne 
o 

E^(y.*)  •  F  (y)e  (16) 

X  X 

where  F^(y)  satisfies  the  equation 


dy2 


+  j^n2(y)  -  sln^ej  F^(y)  -  0  (17) 


The  effective  dielectric  constant  can  be  rearranged  by  first  dividing  Equation 
(Aa)  by  me 


me  m 
0 


L  ti)  J  j  L  0)  J 


The  Index  of  refraction  n(y)  can  now  be  written  as 

2 


Let 


and 


n2(y)  -  ■  1  - 

0 


^  r^; 

-  J  ^  c  L  m  , 


ra  J  ;r 


1  + 


r  T  ^  w  r  1* 

[^]  ‘*[^1 


‘m  (y) 

u(v)  -  * 


[■ 


1  -  J  “ 

m 


1  + 


_1 -  +  j  -liiLj - - 

m  ■•&) 


(18) 


(19) 


then 


n2(y)  *  1  -  au(y) 


(20) 
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Equation  (17)  can  nuw  ba  virltten  aa 


d^F  (y) 

— [cos^e  -  au(y)  ]  F^(y)  -  0 


The  parabolic  electron  distribution  shown  in  figure  2  can  be  represented  by 
the  equation 


whera 


u(y)  -  ^  (y-y^)’’ 

0  .  [jssil 
m  I  “  J 


'0 

y— ► 


Figure  2,  Parabolic  Electron  Density  Distribution 


Now,  let 
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aU 
-k2 


and 


n  ■  y  - 


Equation  (21)  becomes 


d^F  (y) 

-dy2  +  +  Bn2)F^(y) 


0 


(23) 
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SECTION  IV 

SOLUTION  TO  THE  WAVE  EQUATION 


The  teak  of  this  section  Is  to  solve  Equation  (23).  which  can  be  solved 
by  using  a  power  8er'>8  technique  or  by  transforming  the  equation  Into  one 
having  a  known  solution.  In  this  study,  Equation  (23)  will  be  transformed 
Into  anothor  equation  which  has  a  known  solution  and  an  Inverse  transformation 
will  be  used  to  obtain  the  desired  solution.  Using  the  transformation 


F^(y) 


6  -  jn^  ^ 

Equation  (23)  can  be  written  as  (see  Appendix  II) 


d^v(6) 

d{2 


■f  (1/2  -  6) 


d6 


(vT-f  1A) 
4 


v(6)  ■  0 


(24) 


Letting 

2  ^ 


then  Equation  (24)  becomes 


(25) 


Equation  (25)  is  a  confluent  hypergeometric  differential  equation  (reference  5) 
of  the  form 


x2y'»  +  (c  -  x)  -  ay  «  0 


(26) 


The  general  solution  to  Equation  (25)  is 

v(6)  -  K^Vj^(6)  +  K2V2(6)  (27) 
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where 


and 


Vj(«)  -  ^F^(|;  j;«) 


;(«)  •  ^‘^Vld 


The  symbol  is  the  hypergeometric  function  and  is  defined  as 

»  U)^x“ 

^F^(a;c;x)  “  1  +  I  -(^nTT 

n*!  n 


The  terminology  (a)^  and  (c)^  are  factorial  functions  defined  as 


(A)  -  A(Afl)  (A+2)  . . .  (A+n-1) ,  for  n  >,  1 

n 


and 


(A)  -  1  for  A  ft  0 

o 


(A)  -  0  for  A 

o 


where  A  is  any  number. 

By  si';u::»tltuting  Equation  (27)  into  Equation  (16)  and  letting 

32  - 

the  general  solution  of  Equation  (23)  Is 

62n2 


F^(y)  -  e 


[KiF^l(n)  +  K2F^2<") 


(28) 


(29) 


(30) 


The  hypergeometrlc  functions  F^j^(n)  and  F^2(n)  and  Its  derivatives  can  now  be 
written  as 

(|)(en) 

F,i(n)  -  1  +  I  TTT - 


2n 


nil  (1)  n. 
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/iT-l  ,  -  -  ...^j^. 
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f 


rj(n) 


.  (2n)(|)  B2“n^"-‘ 

V  _ 5 - 

n-l  nl 


e2n^2n-l 


(32) 


Fj(n) 


en 


-  \  2  / 

1+  I 

n»l 


^X2(''> 


6+  I 

n-1 


(f)  •' 


(2n+l)(^)  B2“+‘n2“ 

W^’ 


(33) 


(34) 


Substituting  Equation  (30)  into  Equation  (16) ,  the  electric  field  in  the  plasma 
is  given  by 


.lini 


E^(y.*) 


*  e 


[•^X 


jk  zsinB 


(35) 


The  magnetic  fields  can  be  obtained  from  Equations  (14)  and  (15) 

e2n2 


H 


Jwu, 


B^  [V^l(n)  +  V^jCn)] 


-  [Vxi(^^  Vx2<''>] 


-ik  zsinS 
o 


(36) 


^  2 

k  sine  -  -jk  zslne 
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SECTION  V 

TRANSMISSION  AND  REFLECTION  COEFFICIENTS 

For  a  plane  wave  impinging  on  a  plasma  slab  as  in  figure  1,  the  electric 
and  magnetic  field  can  be  described  by  the  following  equations  (reference  T, : 

For 

y  1  0 

-jk  (y  cose  -  z  sine)  jk  (y  cose  +  z  sin0) 

E^(y,z)  -  E^  e  °  e  °  (38) 


k  cose 


-jk  (y  cose  -  2  sine) 
e 


For 


jk  (y  cose  +  2  cose) 
e 


(39) 


y  1  2y 


0 


-jk  (y  cose  -  2  sine) 

E^(y.z)  “  ®  ° 


(40) 


H^(y.z) 


k  cose  jk  (y  cose  -  z  sine) 

_ E  e  ° 

T 


(41) 


The  boundary  conditions  at  y  -  0  can  be  applied  by  equating  Equation  (38)  to 
Equation  (35)  and  Equation  (36)  to  Equation  (39),  and  noting  that  at  y  *  0, 
n  *  y^.  This  will  give 


'  [VxlK)^  Vx2(-yo)  ] 


(42) 
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2 


Ik  cos  8 
o 


+  •'2^x2{-yo)  ] 


(43) 


The.  hypetgeometrlc  functions  K^i(n)  snd  y^J(n)  of  Equations  (31)  and  (34)  are 
even  functions  of  n.  ^^[(n)  and  of  Equations  (32)  and  (33)  ate  odd 

functions  of  n.  Therefore, 

^xl(”^o)  “  ^xl(^o) 

’'xibo)  •  -^x.l(yo) 

^x2(-yo)  ■  -^x2(yo) 

^x2  (-^o)  •  ^x2(yo) 

Equations  (42)  and  (43)  can  now  be  written  as 

+  '  [Vxll^o)  -  Vx2(yo)  ]  (44) 

[Vxi(yo)  -  Vx2(yo)  ] 

-  [•^I’^xK^o)  -  '^2"xHyo)  ] 


2 


*  jk  cose 
o 
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The  boundary  conditions  at  y  "  2y^  can  be  applied  b/  equating  Equation  (35)  to 
Equation  (40)  and  Equation  (36)  to  Equation  (41),  and  noting  that  n  -  y^. 

This  will  give 


-2jk^co89  -  -y 

e  -  e 


[  +  ‘'2"x2(yo)  ] 


-21 k  COS0  2 

•*  o  e _ 

Ik  COS0 
o 


[Vxi(M  Vx2{yo)  ] 


-  [Vxi{yo)+ Vx2(yo)] 


Subtracting  Equations  (46)  and  (47)  gives 


\i{yc) 


jk^cos0  I  jk^cosO 


f  \  ^^xl(yo) 

^2  ^x2{yo)  I  ^  ”  jk^cosS  I  jk^cose 


-  0  (48) 


The  ratio  K2/Kj^  can  be  written  as 


^xi(yo)  y  ■  jk^ 

/ 

^x2(^o)  I  jk  cos9 


jCOS0  I 


Uyp) 

jk  COS0 
jk  COS0 

o 
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‘  v-«r«p^«P .-5..  i!v>» 


Solving  for  E^  find  from  Equations  (42)  and  (43)  gives 


6^1 


K, 


'xl(yo) 


/,  8^0  \  ^i(yo) 

\  JkjjCOse  J  jk^cose 


-  K, 


\  ■  JV“®  /  J'^o 


2(^0) 


COS  6 


8^^ 


-  K, 


*^1 

Vpo)  (  ‘  * 

\ 

"xK^o)  ■ 

Jk^cos'e  1  ~ 

jk  C089 

0 

'l*  'i 

^  +  jk^cose  1 

jk  cos 6 

0 

Adding  Equations  (50)  and  (51)  to  Equation  (48)  gives 


Ei-Ki 


^xi{yo) 


(i--^ ) 

\  jkocose  J 


^xl(yo) 

Ik  CO80 
0 


xli  0)  -  \2{^o)  Jk  cose  ■*■  K,  Jk  cosB 

1  '  '  0  1  o 


Solving  for  E^  in  Equations  (46)  and  (47)  gives 


p» 

>1 

j2k  COS0 

».  *^0 

■  e  < 

"1 

-in  J  *_  %  1  1  1 

\  •'xi(yo) 

^xl(^o)  y  jk^cos0  ^ 

/  "  jk^cos0 

- 

- 

(50) 


(51) 


(52) 


(53) 
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^2  |^^x2(yo)  jk^cose 

Adding  Equation  (54)  to  Equation  (48)  gives 


jk^cosS 


jk  CO80 
o 


i2k  COS0 

r.  o  V 

Efp  *6  K 


K 

1  ^x2(yr  ' 


The  power  reflection  and  transmission  coefficients  are  defined  by,  respectively, 

I  V  i  2  I  17  1 2 


2 

Er 

T  - 

and  R  « 

Dividing  Equations  (54)  and  (53)  by  Equation  (52) ,  the  power  transmission  and 
reflection  coefficients  are  given  by,  respectively. 


\i{yo)  *  ^x2{yo) 

^  ~  jk^cos6  ^  ^xl(^o)  ^  jk^cose  ^xl(^o) 


6^yo  ^2  r  1  1 

(^o)  ”  jk^cos0  j^^x2(^o)  ^  jk^co8  0  ^x2(^o)  j 


CO80  x2l 


^  ”  jk^co80  ^  ^xl(^o)  ^  jk^cos0  ^xl(^o) 
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SECTION  VI 

RESULTS  AND  CCRiCLUSIONS 

The  power  trensmlsslon  and  reflection  coefficients  were  prograoned  for  a 
CDC-6600  Digital  Computer.  The  approximation  of  a  parabolic  electron  density 
distribution  for  a  rectangular  glow  discharge  tube  was  checked  and  the  results 
are  shown  in  figure  3^  The  experimental  points  were  obtained  by  personnel  at 
McDonnell  Aircraft  Corporation  (reference  4).  The  transmission  and  reflection 
coefficients  were  also  calculated  using  a  numerical  technique,  after  the 
electron  density  distribution  was  obtained  from  Langmuir  probe  measurements. 

As  can  be  seen  in  figure  3,  the  parabolic  electron  density  approximation  agrees 
very  closely  with  the  experimental  results  and  the  results  obtained  by  the 
numerical  method. 

The  power  transmission  and  reflection  coefficients  are  plotted  as  a 
function  of  signal  frequency,  collision  frequency,  slab  thickness,  and  angle  of 
incidence  in  figures  4  through  11.  Some  general  conclusions  of  the  electro¬ 
magnetic  transmission  properties  of  a  parabolic  plasma  slab  are 

a.  The  plasma  becomes  more  transparent  to  electromagnetic  signals  as  the 
signal  frequency,  in  radian,  increase  above  the  peak  plasma  frequency  (figure  4); 

b.  For  a  plasma  with  collision  frequencies  lower  than  the  peak  plasma 

frequency,  the  reflection  coefficient  changes  abruptly  when  the  signal  frequency 

in  radian  equals  the  peak  plasma  frequency  (figure  5).  'fherefore,  the  peak 

electron  density  of  a  parabolic  plasma,  with  — j — 2,  can  be  determined 

“^pV^o) 

from  reflectometer  measurements; 

c.  When  the  collision  frequency  is  much  higher  than  the  peak  plasma 
frequency  (greater  than  10),  the  transmission  coefficient  is  independent  of 
the  signal  frequency  (figure  6); 

d.  For  a  given  plasma  thickness,  the  transmission  coefficient  is  greater 
for  higher  signal  frequency  (figure  10); 

e.  The  plasma  becomes  more  opaque  to  the  electromagnetic  signal  as  the 
angle  of  signal  incidence  ir  reases  (figure  11); 
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f.  The  transmission  and  reflection  coefficients  for  the  parabolic  electron 
density  case  behave  in  a  manner  similar  to  the  homogeneous  plasma  case,  when 
the  electron  density,  collision  frequency,  signal  frequency,  thickness  and 
angle  of  incidences  are  varied.  However,  the  reflection  coefficient  of  a 
homogeneous  plasma  is  somewhat  higher  because  of  the  sharp  changes  in  the 
boundaries. 
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Figure  3.  Transmission  and  Reflection  Coefficients  for  Glow  Discharge  Tube 
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Figure  9.  Reflection  Coefficient  versus  v/w 
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-  J  - 

Th«  solution  to  Equation  (I~6)  is 

-  Ke  *  (1-8) 

From  the  geometry  shown  in  figure  1,  ■  -k^sinS, 

If  the  condition  is  made  that  the  wave  propagates  in  the  negative  z-direction, 
the  positive  sign  of  the  exponent  can  be  dropped  (reference  7).  The  solution 
to  Equation  (1-7)  can  now  be  written  as 

jk  zsinS 

G  (z)  -  Ke  °  (1-9) 

X 

Substituting  the  solution  for  G  (z)  of  Equation  (1-9)  into  Equation  (1-2) 

X 

will  give 

jk  zsinB 

E,(y.z)  -  F,(y)e  °  (1-10) 

X  X 
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APPENDIX  II 

TRANSFORMATION  INTO  A  CX)NFLUENT  HYPERGEOMETRIC  EQUATION 


Equation  (23)  Is 


d^F  (y) 

—^r-  +  lA  +  Bn2]  F^(y)  -  0 


(II-l) 


Let 


and 


Fjj(y)  “  e  v(6) 


«  -  Jn^ 


Using  the  chain  rule  for  differentiation 


dF 

_ X  ^  _ X  ^ 

dy  "  d6  dy 


(II-2) 


^  -6/2  fdvi^  _  1 

d6  [da  2 


] 


(II-3) 


(II-4) 


[6>/2  ^  v(6)] 


(II-5) 
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Once  oore  applying  the  chain  rule 


dy*  d«  \  <ly  J  ^ 


(II-6) 


d_ 

d6 


(|i)  .  j  .'/*  (i  -  «■«) 

,  1  ^  j  .-«/=  , 


d^F 
_ r 


4b1/2j 


-6/2 


(II-8) 


Substituting  Equation  (II-8)  and  the  expression  for  F  into  Equation  (23)  and 
noting  that  6n^  ■  6,  gives 


« 


Now,  let 


h  - 


and  Equation  (II-9)  becomes 


1A 

z/b* 


(II-9) 


6d^v(6)  ,  /  I 
d«2  \  2 


^  -  T  x*’ 


0 


(II-IO) 
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